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V.O Abstract 

o 

We study, using a gravity dual, the finite temperature dynamics of SU(N C ) gauge 
theory for large N c , with fundamental quark flavours in a quenched approximation, in the 
presence of a fixed charge under a global current. We observe several notable phenomena. 
There is a first order phase transition where the quark condensate jumps discontinuously 
at finite quark mass, generalizing similar transitions seen at zero charge. We find a non- 
zero condensate at zero quark mass above a critical value of the charge, corresponding to 
an analogue of spontaneous chiral symmetry breaking at finite number density. We find 
that the spectrum of mesons contains the expected associated Goldstone ( "pion" ) degrees 
of freedom with a mass dependence on the quark mass that is consistent with the Gell- 
Mann-Oakes-Renner relation. Our tool in these studies is holography, the string dual 
of the gauge theory being the geometry of N c spinning D3-branes at finite temperature, 
probed by a D7-brane. 
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1 Introduction and Conclusions 



There is great interest in understanding the phase structure of Quantum Chromodynamics 
(QCD). There are various regions of the (/ig,T) phase diagram (T is temperature and fiB is a 
chemical potential for baryon number) which are extremely hard to model theoretically using 
traditional field theory techniques. Some of these regimes are being probed experimentally, and 
so there is additional interest in obtaining better theoretical understanding from an immediate 
phenomenological point of view. 

For almost a decade, progress in understanding the dynamics of extended objects in string 
theory has enabled several new techniques-and a fresh and powerful new perspective-to be 
brought to the study of strongly coupled gauge theory. The techniques were sharpened further 
(for large number of colours, N c ,) with holographic examples such as the AdS/CFT correspon- 
dence, and the study of gauge/gravity duals (or, more properly, gauge/string duals) has proven 
to be a fruitful and often elegant pursuit. 

Not long after the birth of AdS/CFT [H El E], where the finite temperature phase structure of 
M = 4 SU (N c ) (at large N c ) was understood in terms of the thermodynamics of Schwarzschild 
black holes in AdS[2j H], it was recognized [5] that progress could be made in understading the 
properties of gauge theories in the presence of a global current, akin to a chemical potential for 
baryon number or isospin, by studying the physics of charged black holes, such as Riessner- 
Nordstromp), Bj, or more general charged black holes with non-trivial scalars[7], in AdS. 

The resulting (/i, T) or (q, T) phase diagrams (where q is conjugate to /i-the analogue of Baryon 
number) of the J\f = 4 guage theory were found to be rather rich, with e.g., a first order phase 
transition line, and even a second order critical point in the (q,T) plane[5l |6j. This led to the 
hopd^that such studies might lead to insights into the physics of the QCD phase diagram. The 
idea is that studies of the effect of the global U(l) symmetry-even though it is not exactly 
baryon number-might lead to physics in the same universality class as the more realistic gauge 
theories, giving insight into QCD at finite temperature and density. 

A more firm footing for this idea should be obtained by the study of those dynamics in the 
presence of fundamental flavours of quark, a key feature of QCD that distinguishes it from 
the M = 4 gauge theory (which of course only has adjoint matter of very specific types and 
quantities allowed by supersymmetry) . While it is notoriously difficult to study this problem 
as a string dual — so far — it is known [H] how to place fundamental flavours into the theory, in 

1 C. V. Johnson, various seminars, 2000 
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the limit N c ^>Nf. This is a sort of "quenched" limit where the quarks feel the gauge dynamics 
but whose presence does not effect the physics {i.e., the quarks do not back-react). 

While this limit cannot capture a great deal of detailed QCD physics, it is hoped that it may well 
provide valuable clues in the study of the phase structure and will certainly be an instructive 
step along the path to fully understanding the case of fully interacting dynamical quarks. 

The global symmetries in question in the case of the M = 4 theory live in the U (l) 3 subgroup 
of the large 5*0(6) C SU(4) R-symmetry enjoyed by the theory. In the string dual this 
corresponds to that fact that the D3-branes (on whose worldvolume the gauge theory resides) 
have an M. 6 transverse to them, In the decoupling limit that gives the AdS / CFT duality, they 
correspond to a subset of the £0(6) isometries of the explicit S 5 of the AdS$ x S 5 spacetime 
where the dual string theory propagates. 

Fundamental dynamical quarks enter the story as D7-brane probes of the geometry. The M = 4 
supersymmetry is now replaced by M = 2 and the quarks lie in hypermultiplets. The IR 6 now 
breaks into an IR 2 transverse to both the D3-branes and the D7-branes, and an M 4 transverse 
to the D3-branes, but internal to the D7-branes. 

In the probe limit in question, the quarks are the ends of the strings stretching from the D7- 
branes to the D3-branes, their masses being set by the distance between the two sets of branes 
in the M 2 . In the decoupling limit (where the N c D-branes have been replaced by AdSs x S 5 
geometry), the location in the plane transverse to the D7-brane controls the quark mass and a 
condensate of quarks, according to[H]: 

lim Liu) ~ m + — H , (1) 

where m = 2na'm q (the bare quark mass) and — c = (tpip)/(8n 3 a'NfTD 7 ) (the quark conden- 
sate), u is a radial coordinate in the decoupling limit, and Nf = 1 in this paper. 

The condensate actually breaks the U(l) symmetry corresponding to rotations in the M. 2 , that 
acts chirally on the ip fields. We explored the allowed values of the condensate — c as a function 
of quark mass m at finite temperature T (where supersymmetry is broken), and with the system 
charged an amount q under a global U(l) current. We did two cases. The first is the chiral 
current associated with rotations transverse to the D7-branes, and the second is a U(l) that 
mixes transverse and parallel rotations. 

We found a rich structure, upon which we will report in this paper. We found that there are 
generalizations of the first order phase transition already observed (and extensively studied, see 
refs.jlOl HU [El [El [El IE]) at q = for this case. In addition, we found that at finite q a new 
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phase appears where there is a non-zero condensate at zero quark mass m. This is an example 
of spontaneous chiral symmetry breaking in the presence of an analogue of finite fixed baryon 
number or isospin charge. 

We study the meson spectrum of our model and identify the massless "pions" (it is really a slight 
misnomer — see ref.JU]) which play the role of the Goldstone bosons of the chiral symmetry 
breaking. We study the behaviour of the pion mass as a function of quark mass m and find that 
it obeys the Gell-Mann-Oakes-Renner relation[16]. Crucially, we also verify that away from 
the range of q where the symmetry breaking is present, the meson has no massless particle at 
m = 0, confirming our identification. 



2 Spinning D3-branes 

To source subgroups of the U(l) 3 C SO (6) global symmetries corresponding to rotations in 
the transverse M 6 , we consider D3-branes with angular momentum. The supergravity descrip- 
tion/metric for spinning D3-branes is[T7J ITS] : 

ds 2 = h; 1 ' 2 j- (i - A) & + M . d2 j + Hi" | - ^ _ ± + r 2 g Hi (dtf + tfdfi) 



where: 



\i=l / d \i=l 



3 2 

i=i 1 



r% sinh (fa) a 3 r\ 
f 2 

Hi = 1 + 



rl = 47rg s Na' 2 = a' 2 r 4 



. 8 / \ 4 



= ,/1+ '2r 3 y \2r 3 
In order to explore the decoupling limit (a' — > 0), which we are interested in, we define: 

r = a'u , rjj = oc'uh, = otqi , 
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such that: 



4 "4 "2 

-4 1 ^ "TaV > S1Ilh (AO -> "tV > 



H, I + 4 , // 



2 ' 123 a' 2 u 4 A ' 



With these definitions, in the decoupling limit, the metric is: 



ds 2 /a' = a* j- (n 1 n 2 n 3 y 1 fdt 2 + r x du 2 + |^ • d#J 

+A-i{^^ (Rdfa - Adt) 2 + R 2 H t di2 2 1 , (3) 



where we have defined: 

,2 



U = f3 , m d A=^* (4) 



There are two different cases that we consider: 



1. A single charge in the direction: qi = q and q 2 = q% = 0. This case is a useful test case 
since it is the simplest problem to explore and provides valuable insight into the next, 
more complicated case. 

2. Three equal charges: q\ = q 2 = q% = q. This case is of greater interest since the super- 
gravity metric, upon compactification, corresponds to the AdS-Reissner-Nordstrom black 
hole [5] (the previous case gives charged black holes with an extra scalar field excited [7]). 



3 Single Charge 

For the single charge case, we choose qi = q and q 2 = q% = in equation ([3]), and the metric 
becomes: 

ds 2 /a' = A^^-H-'fdf + ^dx-dx^+A^if-V + tfde 2 } 

+A-V2 {ft gin 2 9 (Rdfa- Adtf + R 2 cos 2 edtt 2 } , (5) 

where 

/ = " = ! + ^ = -'«+W S . (6) 
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Although there is a change of variables that allows us to rewrite this metric in a more conven- 
tional form, it will be computationally convenient to work with the metric in equation ([5]). In 
order to introduce dynamical quarks into our gauge theory, we probe the geometry with D7- 
branes [8]. We choose the coordinates (9, 0i) to be the transverse coordinates to the D7-brane. 
However, we will define radial coordinates in the _R 6 -subspace transverse to the D-3 branes via: 

p = u cos 9 , L = u sin 9 . (7) 

Although these radial coordinates are not transverse to each other, the non-zero Gl p component 
of the metric vanishes at the asymptotic boundary. Therefore, by expanding: 

L = m+- + o(-\ , (8) 



u 2 \ ui / 

we can still make sense of the coefficient m as being the asymptotic separation of the D-3 and 
D-7 branes. This will allow us to extract the bare quark mass and the condensate value. It 
should be noted however that, unlike the zero charge case considered in refs.JH.j 13j, we cannot 
directly identify the coefficient c with the quark condensate. The correct identification can be 
made after an appropriate change of variables that restores the S 5 symmetry of the geometry 
at the asymptotic boundary. We will explicitly show this later. 

We choose an ansatz 9 = 9(u) and c^0i = 0. With these choices, the worldvolume of the 
D7-brane is given by: 

1/2 



V 7 - det g = (det S 3 ) 1/2 u 2 cos 3 9 



i- (fq 2 cos 2 9 + {f- A 2 H 2 sin 2 9) u 2 ) (l + fR 2 (9') 2 ) 
J ri 



(9) 

From here, we see that the worldvolume goes to zero for two particular cases: 

1. 9 — 7r/2. This case is akin to the Karch-Katz like solutions described in ref . [1 lj . 

2. | fq 2 cos 2 9 + (/ — A 2 7i 2 sin 2 9) u 2 1 = 0. This is akin to the condensate solutions described 
in ref. [TT]. Note the absolute value sign in @. This means that although the world- 
volume goes to zero at this point, the D7-brane does not actually end there. The vanishing 
is a result of the coordinates we are using; the D7-brane continues through and terminates 
on the event horizon defined by G uu = 0. We are only interested in D7-brane embeddings 
that are smooth (no kinks) across this locus, as they correspond to the lowest energy 
solutions. This condition of smoothness, which has only been checked numerically, is 
only satisfied for a particular boundary condition for our scalar fields on this locus, which 
we describe below. 
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We solve the equation of motion for 9 numerically using a shooting technique. We use R = 1 
and uh = 1 throughout our numerical calculations. Such a choice corresponds to measuring 
distances in units of the AdS radius and to fixing the mass of the black hole, respecively. It 
has become customary in the field [TUl Q2J [15] to shoot from the event horizon towards infinity. 
Using IR boundary conditions for the uncharged case has meant that the correct UV boundary 
values are achieved automatically, and the bare quark mass and condensate are extracted in this 
way. For the charged case, because the ergosphere and the event horizon do not coincide, this 
introduces a new complication. For 9 < tt/2, near the event horizon, the equation of motion 
for the embedding resembles: 



/ {A 1 9"(z) + A 2 9\zf) + -= (A 3 9'(z) + AO 

v J 







(10) 



where A are functions of z = u 2 and 9(z). Since / — » at the event horizon, we use equa- 
tion (10) to obtain the necessary boundary condition, which is depicted in figure [I] However, 
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Figure 1: Boundary condition at the event horizon for q = 0.5. 



if we use this boundary condition to shoot to infinity, we find that 9'{z) diverges at the locus 
of points where the worldvolume vanishes. At first glance, this may not be a problem, since 
it is a necessary condition if we wish to have a regular, non-zero worldvolume across the lo- 
cus. However, in this case, the worldvolume diverges at the pseudo-horizon. This results in 
the D7-brane embedding to curve back on itself, as shown in figure [2j We argue that these 
solutions do not correspond to physical D7-brane embeddings, and, therefore, the boundary 
conditions specified at the event horizon are not the ones we are interested in. We get similar 
results if we try to force the worldvolume to be regular and non-zero across the same locus of 
points. Therefore, we impose a new "physical" condition on our embeddings; we want 9'{z) to 
be at least continuously different iable across the locus of vanishing worldvolume. In order to 
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= u cosQ 



Figure 2: Shooting from the event horizon (q = 0.5) results in embeddings curving back . The 
left-most (black) line is the event horizon; the right-most (red) line is the vanishing worlvolume 
locus. The (blue) curve in between is the embedding solution. 



figure out the appropriate boundary conditions at this locus of points, we study the equation 
of motion near the locus: 



F (B 1 9"(z) + B 2 6\z) + B z e\zf + B A 6\zf) + -^L (B 5 9\z) + B 



-0 



(11) 



where Bi depends on z and 0(z), and C\ is a strictly positive function of z,6\z) 2 . Since 
F — > at the locus of vanishing worldvolume, we use this equation to obtain the necessary 
boundary condition, which is depicted in figure 13} Therefore, for 9 < n/2, we shoot from the 



9 ' 
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Figure 3: Boundary condition at the pseudo- horizon for q = 0.5. 

locus described in point 2 above-which is not where the D7-brane ends, in knowledge that 
we are able to "sew" the rest of the D7-brane that ends on the event horizon. Two examples 
of this sewing process are shown in figure [4] to illustrate that it can be done, and we refrain 
from plotting the sewed solutions on future plots. We use the boundary conditions depicted 
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in figure Km For 6 = 7r/2, this corresponds to shooting from the location described in point 1 




= U COS0 



(a) (b) 

Figure 4: Sewing of solutions for one single charge (q = 0.5) in the transverse direction. The 
black (inner most) curve corresponds to the event horizon, the red (slightly above) curve to the 
shooting locus, and the blue (extends to the right) curve to the D7-brane embedding. 



above, where the D7-brane has ended. The boundary condition is given by: 



7T 



/2 



fl'Lc/a = 00 • ( 12 ) 

These boundary conditions were also used in reference [12] ; it is a result of requiring the absence 

of conifold singularities as the S 3 that the D7-brane wraps shrinks to zero size. 

Figure |5(a)| shows examples of the condensate and Karch-Katz-like solutions as they approach 



their respective end points. It is important to notice that in figure 5(b), the Karch-Katz-like 
solutions have a slope of q 2 /2 for large m, so this choice of the constant c is not proportional 
to the condensate, as we remarked earlier. This is a result of our choice of coordinates in the 
expansion of L as given in equation (J8J), for which: 

q 2 

(iH>) ~-c-—m. (13) 

In order to see this, we recall that for the zero charge case (i.e., the AdS Schwarzschild black 
hole studied in ref.JTTJ), as u —>■ 00, we have L ~ (m + -%); the constant c is proportional to 
the quark condensate as shown by ref. [9j. The key geometric difference between the charged 
and uncharged case is that the S 5 gets distorted due to the presence of the charge q. In order 
to be able to identify the condensate, we want to restore the symmetry of the S 5 as much as 
possible as we move away from the horizon. 
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Figure 5: Solutions for one single charge (q = 0.5) in the transverse direction 

By looking at the distorted S 5 part in equation ([3]), one can see that, as we start moving away 
from the horizon (i.e. when we see the effect of % less and less), it starts conformally looking 
like: 



(14) 



8=1 



where 



/ii = sin 9 , /i2 = cos 9 sin if> , /i 3 = cos 9 cos ip ■ (15) 
Following ref.pjJ], the change of variables that would restore the symmetry is given by: 



r 2 + q\ Hi cos 0i , y 2 



r 2 + q\ [i! sin 0i , y 3 = J r 2 + q\ ji 2 cos 2 , 



Vi = 

Va = \]r 2 + { £ [i% sin <f> 2 , Vs = yj 'r 2 + q% /i 3 cos (f> 3 , y 6 = yj 'r 2 + q 2 /i 3 sin 3 . 
In view of this, for the single charge case, one defines the following radial coordinate: 



w 



+ ^u 2 (u 2 + q 2 



u 



H 



q 2 + 2u 2 + 2^u 2 (u 2 + q 2 ) -11 



(16) 



(17) 



1-1 



This allows us to rewrite: 



f- L du 2 + R 2 d9 2 = PL 1 



dr 7 



r 2 + q^ 



+ d9' 



We can now define new coordinates: 



L = yr 2 + q 2 sin 9 , p = r cos 9 , 
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so that we have: 

R 2 (^% + de 2 ) = R2 (dl 2 + dp 2 ) . (19) 

Then from the asymptotic behavior of r for large u: 

V^ + ? = «+^ + ... , (20) 
we can argue that L has the expansion: 

L = m+ 5-^ + ... • 21 

P 

Using the standard argument from ref.[9] we can show that: 

<V^)~-c-my. (22) 

So, in these coordinates that restore the symmetry of the S 5 , the naive condensate value gets 
a constant subtraction. Henceforth, all future plots of — c will subtract off this term. The 
adjusted solution and the phase transition region are shown in figure [6} 



There are several important things to notice. First, in figure 6(a) , unlike the AdS Schwarzschild 
case, we begin to see the condensate becoming negative for small bare quark mass. In addition, 



we see the appearance of apparently negative masses (see figure 7(a) ) . This of course corre- 
sponds to negative L = usm9 value for p — > oo. This effect is to be interpreted as the constant 
coordinate <f>\ flipping to 7r — <f>\. This has been checked by extending 9 to negative values and 
observing that we get the same physics, continued to the negative mass side. This is simply the 
familiar case of realizing that a negative radial coordinate should be flipped back to positive 
values after taking into account the range of the angular coordinates. Here, by considering the 
fixed 0i ansatz, we are required to extend the definition of 9 to negative values. Therefore, 
these solutions would have physical values given by m — > — m and c — > — c. 

Therefore, we see the appearance of a finite condensate gap between the (m = 0, — c = 0) 
solution and the (m = 0, — c ^ 0). Our equal area law, discussed in ref.[10], tells us that 
the lower branch, corresponding to negative values of the condensate, has the lowest energy. 
Therefore, the system chooses the (m = 0, — c 7^ 0) solution. This means that the geometry 
breaks U(l) chiral symmetry spontaneously. We plot the value of the condensate at zero bare 
quark mass in figure |7(b) We see that there is a minimum amount of charge that is needed 
before chiral symmetry is broken. We also notice that the condensate decreases toward zero 
asymptotically as q — ► 00. 
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(c) Condensate in the transition region (d) Energy in the transition region 
Figure 6: Solutions for one single charge in the transverse direction of the D7-brane for q = 0.5 
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Figure 7: Behavior near m = of the condensate for a single charge. 
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4 Three Equal Charges 



When we set all three charges to be equal[5], (i.e., q\ = q 2 = q% = q in equation O)), the metric 
is given by: 

v 2 

ds 2 /a' = -H- 2 fdt 2 + Hf- l du 2 + H—dx-dx 

R 2 

3 

+ ^2{R 2 d/4 + /4 (Rd<pi - Adtf) , (23) 



i=l 



where 



' 7 R 2 u 2 R 2 ' i?(w 2 + g 2 ) ' 1 J 

From a computational point of view, it is again more convenient to use 

L(u) = usin9 = m+ ^ + ( \) . (25) 

u \u 4 y 

However, just like in the single charge case of last section, the Karch-Katz-like solutions exhibit 
a positive slope of q 2 /2 for the condensate, which is again an artifact of the coordinates. To 
show that, we introduce a new radial coordinate r(u) satisfying: 

Hf^du 2 + R 2 = ^(dr 2 + r 2 d6 2 ) , (26) 
and the corresponding radial coordinates are: 

L = r sin 9 , p = r cos 6 . (27) 



r{u)=u + ^- + o{\ ) , (28) 



Using that r{u) from equation (26) has the asymptotic behavior: 
one can show that: 

L = m+(c + m*P) J 2 +0<^ . (29) 
Note also that one can write: 

Ur X du 2 + R 2 = T f . {dp 2 + dl 2 ) , (30) 
1/ + pr 

and following the standard argument presented in ref.[9] one can show that: 

q 2 

~ -c-m— . (31) 
13 



Therefore, to explore the functional dependence of the condensate on the bare quark mass, we 
again choose 9(z = 1/u 2 ) and <pi to be our transverse coordinates to the D7-brane. With this 
choice, the worldvolume of the D7-brane is given by: 



V-det 
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(det S3 



,1/2 cos 



2z 3 



j (f-A 2 H 2 sin 2 6) (n + AR 2 f, 



1/2 



(32) 



We see that when 9 = n/2, the worldvolume goes to zero. This corresponds to the Karch- 
Katz-like solutions that we saw for the AdS Schwarzschild black hole and our single charge 
case. The other case where the worldvolume goes to zero is when: 



R 2 z (f - A 2 H 2 sin 2 9) = l + z (3q 2 - b A z + 3q 4 z + q 6 z 2 ) - b 4 q 2 z 3 sin 2 9{z) = . 



(33) 



This is a similar situation as for the single charge case, where the worldvolume goes to zero, but 
the D7-brane continues through and ends only at the event horizon determined by G uu = 0. We 
can solve the equation of motion for 9 numerically (using R = uh = 1) using the same boundary 



conditions as in equation ( 12 ) for the Karch-Katz-like solutions and the curve in figure 8 for 



9 < 7r/2, and we illustrate the two types of solutions in figure |9j We see exactly the same 




0.25 0.5 0.75 1 1.25 1.5 



Figure 8: Boundary condition at the pseudo-horizon for q = 0.5. 



behavior as in the single charge; for large m, the slope of the — c vs. m plot, (figure 9(b)) is 
given by q 2 /2, as can be expected from the arguments presented above. Therefore, as in the 
single charge case we will henceforth subtract this q 2 /2 contribution from all future — c values 



displayed in our figures. Figure 10 shows solutions for a variety of charges. We can determine 
where the phase transition between condensate and Karch-Katz-like occurs from the equal- 
area law that we determined in ref. [TU]. This allows us to avoid calculating the energy of the 
D7-brane. We again notice that, beyond a critical value of q, we begin to see solutions with 
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0.5 1 1.5 2 



(a) (b) 

Figure 9: Example of solutions for three equal q = 0.6 charges. The event horizon defined by 
quu = o is shown as the black disk. Again, the red solid lines correspond to Karch-Katz-like 
solutions and the blue dashed lines correspond to the condensate solutions. 



negative condensate value (see figure 10(a)); we saw this feature in the single charge case as 
well. 

Notice also that we again see the appearance of solutions with apparently negative bare quark 
mass; this of course corresponds to negative L = usinO value for p — ► oo, an effect seen in 
figure 9(b) and can be explained as before. For those values, one simply swops m — > —m. 

Therefore, we see the appearance of a discrete condensate gap between the (m = 0, — c = 0) 
solution and the (m = 0,-c ^ 0). Using our equal area law, it would tell us that the lower 
branch, corresponding to negative values of the condensate, has the lowest energy. This means 
that our gauge theory will choose to have a non-zero condensate at zero bare quark mass. 
This means the geometry breaks U(l) chiral symmetry spontaneously. As q approaches its 
maximum value (q pa 0.62) after which there is no event horizon, we begin seeing a multiplicity 
of condensate solutions for small bare quark mass, resulting from the rather beautiful spiral of 
solutions (see the first quadrant in figures [II(aJ1 and [lT(b)j Using our equal law again, we 
learn that the lowest branch of the spiral, corresponding to negative values of the condensate, 
has the lowest energy, and as we go "deeper" into the spiral every successive branch has higher 
energy; The other zero-mass condensates are simply metastable states. Therefore, again at zero 
bare quark mass, we have a non-zero negative condensate value and chiral symmetry breaking. 



We plot the condensate value for zero bare quark mass as a function of the charge q in figure 12 
We see again that there is a minimum q below which chiral symmetry is unbroken. 
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(c) q = 0.62 
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Figure 10: Solutions for three equal charges for various q. The dashed green line represents 
where the phase transition occurs. (In the dramatically magnified regions shown on the right, 
there is some jaggedness of the curves, as we are near the limit of our numerical accuracy. This 
should not be a cause for concern.) 



5 Meson Spectrum 

In order to confirm that we are indeed seeing chiral symmetry breaking, we calculate the 
quadratic fluctuations in the transverse coordinate <p% about the embedding solutions for 9(z) 
we found earlier. In the gauge theory, this will correspond to a pseudo-scalar meson, as described 
in ref . [9] . As a reminder, in ref . [20] , the exact meson spectrum for the AdSs x S 5 background 
was found to be given by: 

2 777 i 

M ( n ^) = T^y / (n + i + l)(n + £ + 2) , (34) 

where I labels the order of the spherical harmonic expansion, and n is a positive integer that 
represents the order of the mode. In our meson spectrum calculations, we have only considered 
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Figure 11: Zoom of the behavior near the origin of the condensate as a function of bare quark 
mass 
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Figure 12: Value of the condensate for the three charge case at zero bare quark mass as a 
function of the magnitude of the charge 



the n = £ = solutions for the fluctuations in 9 and (p. We will choose an ansatz of the form: 

5(pi = x( z ) ex P (—iwt) 1 (35) 

where Re [a;] = M corresponds to the meson's 4-dimensional mass and Im[u;] corresponds to the 
inverse of the meson's lifetime. We then recalculate the DBI action of the D7-brane probe and 
expand the action to second order in xi z )- The classical equations of motion for the fluctuation 
can then be derived using the Euler-Lagrange equation. The first thing we find is that near the 
event horizon, the equation of motion does not have the form of a Heun equation; this means 
that we do not have quasi-normal modes at the event horizon. This means that the condensate 
embeddings and the mesons associated with their fluctuations are stable. This is the opposite of 
what has been found for the AdS Schwarzschild black hole background [2"Tl I2"2l ITU] . Therefore, 
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we only need to consider the real part of u in our analysis. In order to solve for the mass of the 
meson, we again solve our equation of motion by shooting from where the worldvolume is zero 
towards infinity. Our boundary conditions for the condensate embeddings is determined by the 
equation of motion, as before to require that the function x{ z ) is smooth across the vanishing 
worldvolume locus. 



X{z) 



-3 



w.v. — >0 



w.v. — >0 



10 

LU 2 <d(6) 



(36) 



where Q(6) is shown in figure 13 Our boundary conditions for the Karch-Katz-like embeddings 
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Figure 13: Boundary Condition function ©(#) for the single charge case 



is simply given by: 



w.v. — >0 



w.v. — >0 



10- 3 

. 



(37) 



Solutions for the single charge case are shown in figure 14 for a case where we do not see chiral 
symmetry breaking and a case where we do see symmetry breaking. We see that for the case 
where we claimed to see a non-zero condensate at zero bare quark mass (q = 0.7), we definitely 



see a massless meson (see figure 14(b) for the case where we had no condensate at zero bare 



quark mass, we have no massless meson (see figure 14(a)). In addition, we plot the small m 



behavior of the meson in figure [15} If this is indeed the Goldstone boson associated with chiral 
symmetry breaking, it should obey the Gell-Mann-Oakes-Renner relation[16], which states 
that: 

M n = yj-m q {i>i>)lft , (38) 
where m q is the bare quark mass, and / w is the pion decay constant. We can fit our curve nicely 
with a square root function as: 

M = 3.57342 m 1/2 , (39) 
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Figure 14: 0i Meson spectrum for the single charge. The dashed line is the solution for the 
AdSs x S5 case. 



which serves as a nice confirmation of our results. 




0.00005 0.0001 0.00015 0.0002 0.00025 0.0003 0.00035 0.0004 

Figure 15: Small bare quark mass behavior of the single charge 0i meson. The black line is the 
m 1//2 best fit curve 



We find we have very much the same results for the three charge case. The equations of motion 
do not form a Heun equation at the event horizon, so we have stable mesons. The boundary 



conditions are the same as in equations (36) and (37), with figure 16 showing the function 0(0). 
Solutions for the three charge case are shown in figure 17 for a case where we do not see chiral 
symmetry breaking and a case where we do see symmetry breaking. We again find a massless 



meson in the case where we have a condensate at zero bare quark mass (see figure 17(b) . In 
addition, we can fit the curve nicely for small bare quark mass (as shown in figure 18) as: 

M = 3.45336 m 1/2 , (40) 
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Figure 16: Boundary Condition function 0(0) for the three charge case 
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Figure 17: (f>i Meson spectrum for the three charge. The dashed line is the solution for the 
AdSs x S5 case. 
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Figure 18: Small bare quark mass behavior of the three charge <p\ meson. The black line is the 
m 1 / 2 best fit curve 



once again confirming that we have chiral symmetry breaking beyond some critical value of q. 
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